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Abstract 

A partial wave analysis is presented of Crystal Barrel data on pp — )• vr'^vr'^, rjrj and rfrj' 
from 600 to 1940 MeV/c, combined with earlier data on da/dQ and P for pp — t- vt^tt^. 
The following s-channel / = resonances are identified: (i) J^*^ = 5 with mass and 
width (M, r) at (2295 ± 30, 235^11) MeV, (ii) J^'^ = 4++ at (2020 ± 12, 170 ± 15) MeV 
and (2300 ± 25, 270 ± 50) MeV, (iii) ^D^^ J^^ = 3"- at (1960 ± 15, 150 ± 25) MeV and 
(2210 ± 40, 360 ± 55) MeV, and a state at (23001^1], 340 ± 150) MeV, (iv) J^'^ = 2++ 
at (1910 ±30, 260 ±40) MeV, (2020 ±30, 275 ± 35) MeV, (2230 ± 30, 245 ±45) MeV, and 
(2300 ± 35, 290 ± 50) MeV, (v) J^<^ = 1~ at (2005 ± 40, 275 ± 75) MeV, and (2165 ± 40, 
16013,^°) MeV, and (vi) J^^ = 0++ at (2005 ± 30, 305 ± 50) MeV, (2105 ± 15, 200 ± 25) 
MeV, and (2320 ±30, 175 ±45) MeV. In addition, there is a less well defined 6^"'' resonance 
at 2485 ± 40 MeV, with V = 410 ± 90 MeV. For every , almost ah these resonances 
lie on well defined linear trajectories of mass squared v. excitation number. The slope is 
1.10 ± 0.03 Gev^ per excitation. 

The /o(2105) has strong coupling to 7777, but much weaker coupling to Its flavour 

mixing angle between qq and ss is (59 — 71.6)°, i.e. dominant decays to ss. Such decays 
and its strong production in pp interactions strongly suggest exotic character. 



PACS: 13.75CS, 14.20GK, 14.40 
Keywords: mesons, resonances, annihilation 

1 Introduction 



In an accompanying paper [1], data are presented from the Crystal Barrel experiment on pp — )■ 
tt'^tt^, rjrj and rjrj' at nine beam momenta from 600 to 1940 MeV/c; this corresponds to the mass 
range 1960-2410 MeV. The objective of the present paper is to add these new results to earlier 
data on pp — )■ 7r~7r+ and do a combined partial wave analysis. Data on da/dQ were reported 
by Eisenhandler et al. [2] and on polarisation P by Carter et al. [3] from 990 to 2430 MeV/c 
(masses up to 2580 MeV). Further data on both da/dQ and P from 360 to 1550 MeV/c were 
reported by Hasan et al. [4]. The polarisation data play a vital role in the analysis, since they 
separate cleanly states with L = J±l, e.g. ^P2 and ^^2. These data help greatly in consolidating 
the analyses of vr^vr^, rjr] and rjrj' data and help identify resonances. The data of Hasan et al. 
are particularly valuable, since they cover the important mass range down to 1910 MeV. 

iNow at IHEP, Beijinj 100039, China 
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There have been several earher amphtude analyses of the data on pp — )■ n~n~^ [^-H]- We 
find quite good agreement with these papers, though we locate additional features, noteably 
four states. The earlier analyses were restricted to the tttt channel only. One of our primary 
objectives is to include rjrj and rjrj' data, so as to examine the SU(3) character of the fitted 
resonances. 

We also use data on 4"*" states observed in pp riir^Tr^ [12]. They provide further constraints 
on the relative sign and magnitudes of ^'H^ and "^^4 amplitudes. 

The analysis reported here is in terms of s-channel resonances up to J^^ = 6~^~^, plus back- 
grounds which are needed only in the low partial waves 0"'""'" and 1 . There are a number of 
reasons for this approach in terms of resonances. Firstly, we find that the well known /4(2050) 
makes a large and unavoidable contribution. It fits to a mass and width slightly lower than 
averages quoted by the Particle Data Group (PDG) [13]. Because we need to make use of its 
fitted mass, we shall hereafter refer to it as /4(2020). Once this resonance is introduced, it 
defines phases over a considerable mass range. Relative phases with other partial waves are ac- 
curately determined and it becomes impossible to avoid introducing further resonances into all 
other partial waves. This is not surprising, since resonances are predicted in this mass range by 
extrapolation from experience at lower masses. A second related point concerns the polarisation. 
It lies close to +1 over much of the angular range at many momenta and inevitably requires 
large imaginary parts in many partial wave amplitudes. The quantity Pda/ dQ is proportional to 
the imaginary part of the interference between partial waves. Differential cross sections measure 
the real parts of the same interferences. In order to achieve consistency between the principle of 
analyticity and these real and imaginary parts of interferences, we find resonance behaviour in 
all partial waves unavoidable. This is what was found in all earlier analyses except that of Kloet 
and Myrher [11], which focussed on the narrow mass range below a momentum of 1 GeV/c. 



2 Procedures and Formulae 

For completeness, we repeat the standard formulae used in earlier work. The differential cross 
section may be expressed in terms of spin-flip and non-flip helicity amplitudes: 



da/dQ = + (1) 

1 Jmax 

= -Y.{2J+l)fi^Pj{cose) (2) 



- if: ^r^^Ji-PKcose). (3) 

4 j=i JJ(J + 1) 



Here PJ^ are Legendre polynomials. Partial hehcity amplitude f"^ are deflned in terms of angular 
momentum partial waves Tl^j according to: 

V2JTT/^+ = v^Tj_i,j - vOTlTj+i,j (4) 
V2jTlfi_ = ^7TTTj_i,j + VJTj+,,j. (5) 
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Our procedure is to express the Tl^j as sums over resonances, up to 5 for each J value: 



Here arc standard Blatt-Weisskopf centrifugal barrier factors for angular momentum L in 
terms of momentum p in the entrance channel pp and q in the exit meson channel. Explicit 
formulae are given by Chung [14]. They guarantee the correct threshold behaviour in the pp 
channel. The Gi are real coupling constants and 0j are phases for each resonance, arising from 
final state interactions. 

We now comment in general terms on this approach. Firstly, resonances overlap and interact 
to some degree, e.g. via common decay channels. In principle, a K-matrix approach is desirable, 
in order to investigate such interactions. In reality, this is impracticable at present, because 
of the large number of open channels, and almost complete lack of information about many 
of them. A practical point is that the speed of convergence of the fitting procedure is helped 
greatly by using as basis states the eigenstates of the scattering amplitudes; these are T-matrix 
poles. We therefore view the resonances fitted to the data as a parametrisation from which 
further analysis can unfold the physics content. A first step is to find how many resonances are 
required and what their properties are. 

A second concern is that t-channel exchanges also undoubtedly contribute. What effect are 
they likely to have, if any? In each partial wave, they lead to left-hand singularities which arc 
very distant; the left-hand cut opens below s ~ 1 GeV^. For J < 3, there are further resonances 
lying between the mass range discussed here and the left-hand cuts. Any i-channel amplitude 
will acquire the phase dictated by s-channel resonances in each partial wave, via rescattering 
within the present mass range. The effect of any t-channel exchange is therefore to introduce 
some slow s-dependence into the resonance width. A classic example of this is the effect of 
the nucleon pole on the shape of the A(1232) resonance, first treated by the theory of Chew 
and Low [15]. In the present mass range, fitting an s-dependent width to each resonance is 
an impractical luxury, because there are too many unknowns. For example, form factors and 
the opening of new channels will distort the shape of each resonance. For simplicity, we adopt 
constant widths. 

One of our objectives is to examine the SU(3) content of resonances. The simplest approach 
would be to require each resonance to have the same phase 0j for all its decay channels. Secondly, 
amplitudes for decay to r]r] and 7777' are in principle related to those for tttt via the quark content 
of the resonance, and via the well-known composition of 77 and 77' in terms of singlet and octet 
states and the pseudoscalar mixing angle © : 



where cos© ~ 0.8 and sin© ~ 0.6 [13]. Amplitudes for decay oi I — Q qq combinations to 7r°7r°, 
7777 and 7777' may then be written compactly as: 




(6) 



qq> — cos©|77 > -|-sin©|77' > 
ss> — — sin©|77 > -|-cos©|77' >, 



(7) 
(8) 




(9) 



\/A[sin^ 0|?7r7 > -|-cos^ ©17777' > —\/2 cos © sin ©| 7777' >]. 



(10) 
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We have introduced into the latter equation a factor A. This factor has been fitted empirically 
in Ref. [16] to a large variety of data on strange and non-strange final states at high energies. 
Its value is -\/A = 0.8 — 0.9 and wc adopt the central value of 0.85. 

The resonances R we observe in present data may be linear combinations of qq and ss: 

= cos$|gg > +sin$|ss > . (11) 

We use qq as a, shorthand to denote {uu + dd) / -\/2. We expect the observed resonances to be 
dominantly qq in view of their production in pp reactions; however, we would like to test this 
against the observed branching ratios to tttt, 7777 and r]r]' . Amplitudes for the decay of R to these 
three channels are given by the three equations: 

= ^ (12) 
f{r]ri) = ^^(cos'e + ^sin'etan^) (13) 

f{r]r]') = cose sin 0(1 - V2A tan*). (14) 

In our earliest attempts to fit data on t]T] and it immediately became obvious that there 
are major problems in applying these SU(3) relations. The best fit we shall show later has a 
of 10585. If one sets the phase angles (t)i of each resonance to the same value in all three decay 
channels tt^tt", 7777 and r/r/' and takes the resonances to be pure qq, increases dramatically to 
about 45000. The fit is visibly awful. Some deviation from strict SU (3) is obviously required. 
It is necessary to introduce some latitude into phase angles (pi to different channels, or into the 
flavour composition of the resonances, or both. 

To examine where the problem lies, we have tried fitting two extreme scenarios. The first is to 
take all resonances to be pure qq by setting all $ to zero; if the phases 0j of all resonances are left 
free, — )■ 14558. This is not a huge increase in over our best value of 10585. However, the 
fit to both rjT] and tt^tt" angular distributions is visibly poor, particularly for r]r]. Furthermore, 
some of the phases 0j depart from zero by unreasonably large amounts > 90°. So it seems that 
some departure from strict SU(3) is unavoidable. 

The other extreme is to set all 0j to zero and attempt to fit purely by adjusting fiavour 
mixing angles $. This turns out to be more successful, giving a of 11527. We shall give 
detailed results below. Most of the mixing angles optimise in the range 0-20°, i.e. close to pure 
qq states. These small mixing angles allow modest changes in branching ratios between 7r°7r°, 
•qr] and r]r]' . However, one resonance, /o(2105), requires considerably larger mixing, with $ in 
the range 59-71.6°. The integrated cross sections shown below for vr^vr", 777/ and 7777' are fairly 
close to the values predicted by simple SU(3) on average. The question therefore is whether it 
is wise to set phase angles 0j strictly at zero. We now consider several considerations bearing 
on this point. 

Firstly, each resonance rides on the tails of other resonances. We have had experience of 
fitting data at lower masses by both T-matrix techniques [17] and the K-matrix [18]. From this 
experience, we have learned that mixing between states commonly gives rise to deviations of 
phase angles (pi betweeen different decay channels in the range -15 to +15°. Even if K-matrix 
poles have the same phases, T-matrix poles can have different phases; the difference depends 
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on the separation of poles, resonances widths and the mixing through decay channels. In view 
of this experience, we feel it unwise to demand strictly the same phases (pi in all channels. For 
most resonances, wc have allowed 0.; to optimise freely for each decay channel with differences 
limited to the range -15 to +15°. In isolated cases, where they seem to need it, this range has 
been allowed to extend to ±30°. 

Angles $ for flavor mixing mostly optimise naturally in the range -15 to +15 °. However, 
for large spin, it is likely that differences in centrifugal barriers and form factors between tttt, 
rjrj and rjrj' may affect branching ratios. Resonances arc commonly believed to have radii of 
0.8-1.0 fm. In tttt decays, the wavelength of each outgoing tt is 1.1 fm for a resonance mass 
of 2.2 GeV. Overlap of wave functions will play a strong role in determining matrix elements. 
Even quite small momentum differences between tttt and 7777 and larger differences for r]r]' may 
lead to significant departures from SU(3). We observe a problem in fitting 7777 cross sections 
simultaneously with vr^Tr*^ unless some departure from strict SU(3) is allowed. Our final com- 
promise is therefore to allow fiavour mixing angles to vary in the range to ±30°. All optimise 
naturally in this range, except for /o(2105), which must be allowed complete freedom. As far 
as resonance masses and widths are concerned, the freedom in fiavour mixing angles fortunately 
has little effect; masses and widths are determined by the locations of singularities, which are 
independent of mixing and corresponding branching ratios. 

A further unknown is the radius to be given to centrifugal barriers. The masses fitted to 4+ 
states are sensitive to this radius, since the centrifugal barrier in the pp channel is strong. If 
the /4(2020) is to be fitted to the PDG value of 2044 MeV, the radius required is unreasonably 
large, namely 1.5 fm. The fitted value is 0.88 fm. 

Ratios of amplitudes for L = J± 1, i.e. rj = are fitted to real constants for 

each resonance. Since this linear combination refers to mixing within the pp channel, it is taken 
to be the same for decays to tttt, 7777 and 7777'. For the /4(2020), one expects the channel with 
L = J — 1 to dominate over L — J+1, because of the centrifugal barrier. This is precisely what 
we find. Likewise, the 5~ state is dominantly ^G^. This mixing is determined very precisely by 
the polarisation data. For 6'^, we assume pure ^Hq. 

The data on 7r°7r°, 7777 and 7777' stop at 1940 MeV/c. Polarisation data stop at 2200 MeV/c. 
From n~7r'^ data above 2 GeV/c, it is clear that a 6"^ resonance is required at 2485 ± 40 MeV. 
A resonance close to this mass has been observed decaying to tttt by the GAMS collaboration 
[19], at 2510 MeV. It seems likely that there will be further high spin resonances of similar 
mass. They lie so close to the top of the available mass range that we cannot establish their 
parameters precisely. Nonetheless, it is convenient to use 5~, 4+, and 2+ resonances in the range 
2500-2620 MeV to parametrise the data at the highest momenta. The tail of the 6+ resonance in 
the mass range LEAR below 2410 MeV needs to be described accurately, since its interferences 
with states of lower spin are important for fitting Legendre polynomials up to order 10 required 
by the tttt and 7777 data. 

3 Results 

We summarise first the expected states. The /4(2020) is well known and there is plenty of 
evidence, summarised by the Particle Data Group, for its radial excitation /4(2300). There 
is a known ps at 2330 MeV. Consequently, resonances having spins with all below 4+ are 
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expected around 2020 MeV and with all below 5~ around 2300 MeV. In the mass range 
1910-2410 MeV, our main focus of attention, this implies two ^D^ qq states plus one with 
near 2300 MeV, two and two states, four 1 states and two 0^. We are able to identify 
all of these except for two 1^ states. It is hardly surprising that it is difficult to disentangle two 
states of low . 

3.1 Uniqueness 

The ambiguities we have encountered concern signs of the amphtude ratios =^ i?4/^F4 and 
r2 F2/^P2- These are the ratios of amplitudes after centrifugal barrier factors are factored 
out, i.e. they are asymptotic ratios as s ^ 00. We have located two solutions with similar x^, 
one with r4 positive and the other with negative. They both contain the same resonances with 
similar masses and widths, and differ mostly in coupling constants. However, analysis of t^tt^tt^ 
data [12,20] allows a clear distinction between these alternatives. Clebsch-Gordan coefficients 
in equns. (4) and (5) are such that positive r4 requires dominance of M = 1 amplitudes for 
74(2300). [Here Af is the projection of the spin of the initial state along the beam direction; 
M = 1 corresponds to in the helicity basis]. This agrees with the t^tt^tt^ analysis. Negative 
r4 is rejected strongly by that analysis. 

There is an argument concerning matrix elements for decay which supports this conclusion. 
For this high spin, the wave function of the resonance is peaked strongly at a radius of 0.8-1.0 
fm. Wave functions for the pp channel are described by spherical Bessel functions. For momenta 
in the mass range under discussion, the first zero of jiikr) lies outside 2 fm for both L = 3 and 
L = 5. Overlap of wave functions with the resonance leads to the expectation that pp ^F^ and 
^H^ will couple to the resonance with the same sign, hence r4 positive. We find that the fitted 
amphtude ratios hl^G^ and Gs/^D^ are likewise positive, as one would expect from 
the same argument. 

We shall present evidence for four 2+ resonances at 1910, 2020, 2230 and 2300 MeV. We 
shall argue for several reasons that the /2(2020) is the first qq '^F2 state. We have examined 
all alternative sign combinations of r2 for these four resonances. The 7r~7r~^ polarisation data 
are sensitive to these signs. The best solution is obtained with signs respectively +, +, -, - in 
the order /2(1910) to /2(2300). The second l)est solution has signs -, -, -, -, but a of 13850, 
i.e. worse by 3265 than the best solution. The increase in is mostly for polarisation data, 
where there is a clear systematic discrepancy between data and fit for the second alternative. 
For 72(2020), r2 is again expected to be positive. The first zero of the Bessel function ji{kr) 
(coupling to ^P2) again lies outside 2 fm; so the pp momentum is low enough that the argument 
given above remains unchanged. 

For /2(2300), which we shall interpret as the radial excitation of /2(2020), the momentum in 
the pp channel is such that the first zero of ji{kr) lies at 1.28 fm. When one allows for (i) the 
fact that the wave function will be attracted to smaller r by the resonant interaction (mixing in 
a component due to the irregular function ni(/cr)), and (ii) the node in the radial wave function 
of the resonance, it is possible to arrive at a negative r2, as observed. For qq ^P2 states, it is not 
possible to make reliable predictions, because of the large number of nodes in the radial wave 
functions of the resonances. They are, however, found to show the same general variation with 
s as for other J^, i.e. positive for low p momenta and going negative for high momenta. 
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3.2 Quality of the fit 



A solution is obtained in typically 1-2 minutes of computing; consequently several thousand fits 
have been made, exploring the systematic effects of including or omitting various resonances. In 
the course of this work, we have found, from the observed variations from fit to fit, that errors on 
resonance masses and widths may be assigned from changes of ~ 20. This is the quantitative 
criterion we adopt in assessing errors given here. It is larger than statistics (Ax^ = 1) and 
therefore conservative. It covers almost all of the variations which have been observed between 
fits with differing ingredients; where it does not, errors have been increased to cover the observed 
variations. 

The Ti^-K^ data have statistical errors which arc much smaller than those of other data. Con- 
sequently, there is the danger that delicate features of the 7r°7r° data cause problems in fitting 
other channels. To remedy this, we reduce the weight of the tt^tt*^ data by a factor 3, and in- 
crease the weight of -qr]' by a factor 2. This does not introduce any major qualitative changes, 
but speeds convergence greatly. 

Values of for all data sets are given in Table 1. They are somewhat above 1 per point, 
particularly for tt'^tt'', because of the weighting mentioned above. It may refiect the presence of 
systematic errors in some data points, but also may reflect small missing components in the fit 
(e.g. backgrounds from the tails of lower mass resonances), which cannot presently be identified 
unambiguously. Our objective is to locate the essential features of the data, and these are stable 
against small variants in the parametrisation of the amplitudes. 



Data 


Author 


Points 


A" 


7r~7r"'" 


Hasan 


1000 


2677 


P 


Hasan 


1000 


1812 


7r~7r+ 


Eisenhandler 


960 


1702 


P 


Eisenhandler 


1078 


1720 


7r°7r° 


Anisovich 


360 


1853 


7777 


Anisovich 


360 


667 


rjr]' 


Anisovich 


144 


142 


Normalisations 


Anisovich 


27 


12 



Table 1: Contributions to of the fit. 



3.3 Normalisation 

Each of the data sets on tt^tt^, rji] and rir]' is given a normalisation constant which has been de- 
termined in Ref . [1] . In our final fits, the normalisation is optimised by allowing each data set to 
vary in normalisation and including into x^ a contribution for the deviation of the normalisation 
from its experimental value. 

As described in [1], the normalisation of Crystal Barrel data for Tr^vr" is on average slightly 
more than a factor 2 higher than that of data of Dulude et al. [21]. This calls into question 
the absolute normalisation of all the neutral data. In order to test this, we have tried scaling 
the normalisation of all Crystal Barrel data for tt'^tt^, rjr] and rjr)' down by this factor 2, so as 



7 



to agree on average with Dulude et al. It results in a huge increase in from 10585 to 16893. 
The polarisation data are particularly sensitive to this change, increasing from 3532 to 7940. 

The reason for this is straightforward, but important. Partial waves in the n^n^ channel have 
isospin and even spin. These amplitudes contributes also to tt^tt^ data, but interfere with 1 = 1 
components with odd spin. This interference gives rise to forward-backward asymmetries in the 
differential cross section, depending on real parts of interferences. Values of Pda/dfl contain 
interferences depending on imaginary parts of identical interference terms. The normalisation of 
Dulude et al. requires much larger 1=1 amplitudes than our data. A fit to the data of Dulude 
et al. alters the phases of interferences, in order to fit asymmetries in differential cross sections. 
But these changes are inconsistent with the polarisation data. This inconsistency rules out the 
normalisation of Dulude et al. If we float the normalisation of Crystal Barrel data freely, it 
optimises at 0.989 times the published values. An increase of of 20 (the criterion we have 
adopted above in discussing resonance masses and widths) corresponds to a change of ±0.023 in 
normalisation. This lies within the normalisation errors quoted for Crystal Barrel data, namely 
±0.03 for momenta of 1050 MeV/c upwards and ±0.06 at 600 and 900 MeV/c. 

In the final analysis, we do not use the factor 0.989. The normalisation of Crystal Barrel data 
is used with its experimental normalisation and with the errors (3-6%) obtained experimentally. 

We include the data of Dulude et al. for the determination of the shape of the angular 
distribution. They provide some valuable points close to cos 9=1. However, their normalisation 
is fitted freely. 

3.4 Comparison of data and fit 

Figs. 1-5 show comparisons of our final fit with data for n~7i~^ differential cross sections, 
polarisation and integrated cross sections. The 7r°7r° data, rjr] and rjr)' are available only up to 
centre of mass scattering angle 6 given by cos^ < 0.875. Their integrated cross sections are 
therefore compared in Fig. 5 with the fit only over this angular range. The 7r~7r"'" data are 
available over the full angular range and are integrated over this range. 

Panels on Figs. 1-3 labelled H refer to Hasan et al; those labelled E refer to Eisenhandler et 
al. or the same experiment of Carter et al. On Fig. 1, the angular distributions are mostly fitted 
well. There are small discrepancies near cos^ = -1 for momenta of 790, 988, 1089, 1140, and 1291 
MeV/c. However, in all cases, data from the other experiment at nearby momenta show either 
no discrepancy or small ones. We therefore take these discrepancies to be due to experimental 
error or statistics. There is some discrepancy with da/dfl near cos 6' = -1 at 467,497 and 585 
MeV. We can find no solution which resolves this discrepancy without rapid changes with mass 
in high partial waves J = 4, 5 or 6; at these momenta, contributions from J = 5 and 6 should 
be small. 

On Fig. 2 there arc some discrepancies at individual momenta, but nothing systematic over 
a range covered by a resonance. Above 2 GeV/c, we are fitting only with tails of resonances 
from lower masses, plus resonances with = 6"*", 5~, 4+, and 2+ in the range 2500-2600 
MeV. It seems likely that there will be further low spin resonances in this mass range. We 
have tried adding one by one resonances with = 3~, 1~ and 0+, but none gives a significant 
improvement in x^- Fo^^ polarisation data on Fig. 3, the agreement between data and fit is 
generally satisfactory and free of systematic trends. 

The fit to 7r°7r° data is shown in Fig. 4 and the accompanying paper on data analysis [1]. It 
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Figure 1: Differential cross sections from 360 to 1300 MeV/c, compared witli tlie fit (full curves); 
panels labelled H are data of Hasan ct al. and those labelled E are from Eisenhandler et al. 
Each distribution is normalised so as to integrate to 27r. 
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Figure 2: As Fig. 1, 1351 to 2230 MeV/c. 
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Figure 3: Polarisation data compared with the fit (full curves). 
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gure 4: Differential cross sections for 7r°7r°, rjrj and rjrj' compared with the fit (full curves), 
these are the average of 47 and 87 data reported in Ref. [1]. 
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Figure 5: Integrated cross sections compared with the fit (full curves). For 7r°7r°, rjr] and rjrj', the 
integration is from cos ^ = to 0.875 and is over only one hemisphere; the rjr]' data are averaged 
over 47 and 87 data reported in Ref. 1. For tt"?!"*" the integration is over the full angular range. 
In (a), black circles are data of Hasan et al. and open squares are those of Eisenhandler et al. 
The dashed curve in (a) shows the contribution from the O"*" intensity. 



is good at most momenta. At 1200 MeV/c, one point at cos^ = 0.875 cannot be fitted under 
any circumstances. The acceptance of the detector falls rapidly around cos^ = 0.875, so we 
suspect this point is affected by some systematic error. The data at 1200 MeV/c were actually 
taken in a run separate from all other momenta, possibly accounting for a discrepancy between 
this momentum and others. 

For 7)7) data, shown on Fig. 4 and in the accompanying paper, the fit is generally good, though 
there are some discrepancies near cos^^ = 1 at 900 and 1940 McV/c. Fits to these regions may be 
improved by relaxing the constraints on phase angles 0j. At 600 McV/c, statistics are low and 
it is possible to obtain some variety of fits, varying the strengths of flavour mixing for /4(2020), 
/2(1910) and /2(2020). Near cos^ = 1, there is a dip in rjr] differential cross sections at 900 
MeV, but a peak at 600 MeV/c. Our fits arc a compromise between these rapidly changing 
features. The fit to either momentum may be improved, but not both simultaneously. For rjr)', 
fits are satisfactory, but statistics are low. 

Integrated cross sections, shown on Fig. 5, are described well. For 7r°7r° in Fig. 5(b), the 
peak at 900 MeV/c is mostly due to /4(2020). For 7777 in Fig. 5(c), the striking peak at 1200 
MeV/c comes largely from /o(2105) and /4(2020). 

3.5 Fitted Resonances 

Fitted resonances are shown in Table 2. We shall comment in detail on most of them. All 
those in the mass range 1900-2400 MeV are established securely, with the exception of the 3 
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Mass M 


Width r 


Ax' 


Hasan M 


Hasan F 




(Mev) 


(MeV) 




(MeV) 


(MeV) 


6+ 


2485 ± 40 


410 ± 90 


1776 


- 


- 


5- 


~ 2500 


^ 470 


112 


2881 


310 


5- 


2295 ± 30 


2351^^ 


2534 


2303 


169 


4+ 


~ 2500 


~ 400 


1305 


2813 


257 


4+ 


2300 ± 25 


270 ± 50 


2549 


2314 


278 


4+ 


2020 ± 12 


170 ± 15 


22382 


(2049) 


(203) 


3- 


2300^^1] 


340 ± 150 


183 






3 


2210 ± 40 


360 ± 55 


368 


2232 


220 


3- 


1960 ± 15 


150 ±25 


2957 


2007 


287 


2+ 


~ 2620 


~ 430 


776 






2+ 


2300 ± 35 


290 ± 50 


2879 


2517 


264 


2+ 


2230 ± 30 


245 ± 45 


2290 


2226 


226 


2+ 


2020 ± 30 


275 ± 35 


2980 


1996 


134 


2+ 


1910 ± 30 


260 ± 40 


2286 






1- 


2165 ±40 


160t^^° 


450 


2191 


296 


1- 


2005 ± 40 


275 ± 75 


1341 


1988 


244 


1- 


(1700) 


(180) 


8444 


(1690) 


246 


0+ 


2320 ± 30 


175 ± 45 


1257 


2321 


223 


0+ 


2105 ± 15 


200 ± 25 


4030 


2122 


273 


0+ 


2005 ± 30 


305 ± 50 


370 






0+ 


(1700) 


1000 


2844 


1745 


238 



Table 2: Resonances fitted to the data; errors cover the full range of systematic variations 
observed in a large variety of fits with varying ingredients. Values in parentheses either lie 
outside the range of masses fitted and describe background amplitudes or are fixed. The fourth 
column gives changes in %' when each resonance is removed from the fit and all other parameters 
are re-optimised. The final two columns compare with results of Hasan and Bugg, Ref. [9]. 
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resonance at 2300 MeV, which is rather weak. The fourth column of the table shows changes 
in when each resonance is removed from the fit and all other parameters are re-optimised; 
this measures the significance of each contribution in fitting the data. Most of these changes arc 
enormous and leave no possible doubt about the presence of these resonances. Our impression, 
from the consistency between alternative fits with differing ingredients, is that any resonance 
contributing a change in > 100 is definitely present. However, below A^^ = 300, the mass 
and particularly the width are hard to establish with confidence. Those resonances contributing 
1000 to the improvement in are very secure in all parameters. 

Argand diagrams are shown for tttt amplitudes on Fig. 6 and those for rjrj and rirj' on Fig. 7. 
Resonance loops are observed in all partial waves. We shall comment in detail later. Intensities 
of all partial waves are displayed in Fig. 8. Dashed contributions are for L = J — 1 and dotted 
for L = J + 1. [Interferences disappear from their sum]. 

There is a strong 4"*" resonance at 2020 MeV and a somewhat weaker, but definite, one at 2300 
MeV; there is tentative evidence for another around 2500 MeV. A strong 5~ state is observed at 
2295 and there is tentative evidence for a higher one at ~ 2500 MeV. There is definite evidence 
for the expected G"*" resonance at 2485 ± 40 MeV. For = 3~, there are two well defined ^D^ 
states at 1960 and 2210 McV. There is tentative evidence for a further ^G^ state at 2300 MeV, 
but it is weak; despite this, its angular dependence is distinctive and it is hard to avoid the 
conclusion that it is present, even if its mass and width cannot be determined well. 

An important result is that there are definitely four 2+ states at 1910, 2020, 2230 and 2300 
MeV; these are expected from the quark model for and ^¥2- For — 1~, four states are 
again expected. There is evidence that two are required, and there are tentative indications 
of a third; however, it is not possible to establish the mass and width of this third state. For 
= 0"*", at least two resonances are definitely required. The one at 2105 MeV appears very 
strongly, particularly in the rir] data, and its mass and width are determined precisely. There is 
a further definite 0+ resonance at 2320 MeV. As we shall comment below, there is probably a 
third O"*" state at 2005 MeV, but its mass and width are not well determined because of overlap 
and interference with /o(2105). We now discuss each resonance in turn, starting with the high 
mass region. 

3.6 The mass range 2500-2600 MeV 

Our main objective is to study the momentum range up to 1940 McV/c, where the extensive 
Crystal Barrel data are available. However, we also fit the momentum range above 2 GeV/c with 
the objective of obtaining a reliable determination of the the tails of the high mass resonances 
below 2 GeV/c. The 7r~7r+ differential cross section data extend to 2430 MeV/c and polarisation 
data to 2200 MeV/c. 

We find definite evidence for a resonance at 2485 MeV, in close agreement with the 
result of the GAMS collaboration [19]. Its contribution to is large, as one sees from Table 
2. We find it desirable to introduce further high mass contributions for 5~, 4+, and 2+, as 
one sees from the changes they introduce into x^, listed in Table 2. However, only the lower 
half of each resonance is within the available mass range. The result is, of course, a strong 
correlation between fitted mass and width, so we are unable to assign accurate parameters to 
these resonances. Furthermore, masses of the high spin states are again sensitive to the radius 
of the centrifugal barrier. 
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Figure 6: Argand diagrams for tt+tt^ partial wave amplitudes fj{s). Crosses mark masses at 
100 MeV intervals beginning at 1900 MeV. The scale is such that the integrated partial wave 
cross section in given in /ib by 27r|/j(s)pp(s), where p(s) is the phase space q/y/s. 
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Figure 7: Argand diagrams for 7777 and 7777' amplitudes. Crosses mark masses at 100 MeV intervals 
beginning at 1900 MeV. The scale is as described for Fig. 6. 
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Figure 8: Intensities from individual for channels 7r°7r°, 7777 and 7777'. Dashed curves show 
contributions with J — J — 1 and dotted curves those with L — J + 1. 
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We now comment on individual resonances and illustrate on Fig. 9 the deterioration in the 
fit when they are dropped from the fit. 

3.7 5- and 4+ 

The 5~ resonance at 2295 MeV is very clear in polarisation data. Figs. 9 (a)-(c) illustrate 
the deterioration in fits when it is removed. The width of p5(2295) is 235I40 MeV, somewhat 
less than the 400 MeV quoted by the GAMS collaboration [22] . However, an inspection of the 
GAMS data reveals the possibilty that they have fitted resonances at 2295 and 2500 MeV by a 
single resonance. On Fig. 6 one sees that the 5~ state at 2295 MeV is dominantly ^G^ and the 
well known /4(2020) is almost purely ^F4. The centrifugal barriers for pp will suppress higher L 
strongly. 

The phase advance of 360° for '^^4 on Fig. 6 requires the presence of two 4+ resonances. A 
second one at 2300 MeV is definitely required, but is fairly weak in the tttt channel, leading to 
a large uncertainty on its width. It is more clearly visible in rjr] data. Earlier fits to tttt data 
[5-10] all gave masses in the range 2300-2340 MeV. Crystal Barrel data on the rin^n^ final state 
show a strong peak in f2{1270)r] at 2320 ± 30 MeV with F = 220 ± 30 MeV [12]. Also recent 
VES data display conspicuous 4+ peaks in cucu at 2325 ± 15 McV with F = 235 ± 40 MeV [23] 
and in 7]%% data with M = 2330 ± 10{stat) ± 20{syst) MeV with F = 235 ± 20 ± 40 MeV [24]. 

3.8 3 states 

The 3~ state at 1960 MeV is strong and very secure. Statistically, the error on the mass is ±7 
MeV and on the width is ±14 MeV. In all fits, it never moves outside the mass range 1950-1980 
MeV. We increase errors to cover all systematic variations. An interesting point is that this 
resonance is definitely lower in mass than the corresponding 4+ states. This demonstrates that 
all resonances in this region are not degenerate in mass. Again, the phase advance of 360° for 
^D^ on Fig. 6 requires the presence of at least two resonances. The higher one agrees well in 
mass with the determination of Hasan and Bugg, namely M — 2232 MeV. Both resonances are 
coupled dominantly to pp ^D^. 

In order to demonstrate that the data really demand the presence of each of these resonances, 
we have dropped them one by one from the fit and examined the discrepancies which emerge 
between data and fit. Figs. 9 (d) and (e) illustrate some of the defects when the lower 3~ 
resonance is dropped; Figs. 9(f) and (g) show defects in the fit when the upper 3~ resonance is 
removed. In both cases, there is a tendency for the remaining resonance to become broad and 
span the mass range. 

A further qq ^G^ state is expected in the mass range around 2300 MeV, close to ^^5. There 
is evidence for its presence as p3(2300) of Table 2. However, it improves only by 183, with 
the consequence that its mass and particularly its width are not well determined. It appears 
almost purely in pp ^G^. The angular dependence of this state is distinct from ^D^, making its 
presence very likely, despite the fact that its parameters are not well determined. This angular 
component requires a mass ~ 100 McV above that of p3(2210). 

We have also tried including p3(1690) as a background amplitude. It has only a very marginal 
effect, so we omit it. The reason may well be that it is inhibited near the pp threshold by the 
L — 2 centrifugal barrier. 
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3.9 2+ states 



An important outcome of the present analysis is strong evidence for four 2+ states in this mass 
range, as expected from qq and ^^2. Earher, there has been evidence from VES and GAMS 
groups for an /2(1920) of rather narrow width, F ~ 90 MeV [25,26]. Recent VES data on uu 
give an improved mass determination of 1937 ± 12 MeV and width F = 150 ± 17 MeV [23]. Data 
on r/7r°7r° from Crystal Barrel [12,20] require 2+ states at (a) M= 2020 ± 50 MeV, F = 200 d= 70 
MeV, decaying dominantly to /2(1270)ry, (b) M = 2240 ± 40 MeV with F = 170 ± 50 MeV, 
appearing as a second peak in /2(1270)r7, and (c) M = 2370 ± 50 MeV with F = 320 ± 50 MeV, 
appearing as a strong peak in decays to a2(1320)7r. 

The inclusion of polarisation data is important in clarifying the components required in the 
fit. The /2(1920) appears very strongly with M = 1910 ±30 MeV, F = 260 ±40 MeV. However, 
it is at the bottom end of the range of data we analyse, and we observe some tendency for 
its mass to drift downwards. We view the determination of the mass by VES and GAMS as 
more secure, but our larger width appears reliable. If this resonance is dropped from the fit, 
increases by a definitive amount, namely 2290. There are large effects on the fit at many of the 
low momenta; effects at two momenta are illustrated in Figs. 9(j) and (k). 

The next 2+ state optimises at 2020 ± 30 MeV with F = 275 ± 35 MeV. In earher fits to r]mT 
data [12] there is the possibility of cross-talk between the two nearby resonances at 1910 and 
2020 MeV, both lying at the bottom end of the available mass range in the neutral data. So the 
confirmation of /2(2020) is important. We find a very large change in x^, namely 2980, when 
/2(2020) is removed from the fit and others are re-optimised. Figs. 9(1) and (m) illustrate the 
deterioration in the fit to some data. 

The present analysis confirms the presence of two further 2+ states at 2230 MeV and 2300 
MeV. When the resonance at 2230 MeV is removed and others are re-optimised, g^ts worse 
by 2290, a very large amount. The effect of removing it is to produce a fit with a single resonance 
at 2280 MeV. The effect of removing the highest 2+ state at 2300 MeV is 2879 in x^- Again the 
effect is to move the 2230 MeV resonance close to 2280 MeV and produce the same bad fit as is 
illustrated in Figs 9(n) and (o). 

The highest of these /2 states, now at 2300 ± 30 MeV, was observed in Refs. [12] and [20] 
at 2370 ± 50 MeV. This shift in mass arises from a subtle interference with /2(2230), revealed 
by the tt'tt"*" polarisation data. For the same reason, the width of the /2(2230) is now 245 ± 45 
MeV, rather than the 170 ± 50 MeV reported earlier. 

All four states appear to be dominantly qq. We shall discuss below the flavour mixing angles. 

3.10 1- states 

The fit is poor without two 1~ resonances. The phase advance for ^Di on Fig. 6 is evidence 
for the presence of two states. However, four are to be expected in this mass range, originating 
from pp ^Si and ^Di. We cannot identify more than two with confidence. 

Figs. 9(h) and (i) illustrate the effect of dropping the 1^ state at 2000 MeV and re-optimising 
the rest. There are significant visible discrepancies in fits to both diffential cross sections and 
polarisation at low momenta. The change in is 1341. It couples strongly to pp ^Di and 
this makes it distinctive in polarisation data in the low momentum range where ^Si dominates. 
The 1~ state at 2165 MeV has a less significant effect, improving by only 450. It couples 
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Figure 9: Illustrating the deterioration in the fits to data when various resonances are dropped 
from the fit: (a)-(c) P5(2295), (d)-(e) P3(1960), (f)-(g) P3(2210), (h)-(i) pi(2005), (j)-(k) 
/2(1910), (l)-(m) /2(2020), (n)-(o) /2(2230) or /2(2300), (p)-(r) /o(2105) and (s) /o(2320). 
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dominantly to pp ^Si. When it is removed from the fit, changes are just discernible by eye. 
There are, however, reports of a f ~ resonance at 2150 MeV from GAMS data on the utt channel 
[22]. This mass corresponds closely to what we observe. The width given by GAMS is large, 
320 ± 70 MeV. The width we observe is smaller, but has a considerable error, 160^70° MeV, so 
there is no discrepancy. 

We find it is also essential to include a strong 1^ contribution peaking at or below the pp 
threshold. It fits well as a contribution from p(1700). It is possible that some of this effect 
originates from the ^5*1 threshold in pp, but a fit using only this threshold is somewhat poorer. 

We have searched for further 1~ resonances at higher mass. There is some improvement in x^, 
with an additional state around 2270 MeV. However, the analysis will not support the presence 
of three 1~ states. 

3.11 0+ states 

One of the very striking features of the 7777 data is a strong peak in the integrated cross section 
of Fig. 5(c) at 2100 MeV. It is not fitted by /4(2020) and /2(2020), whatever the centrifgual 
barrier radius. These states with non-zero spin are separated by their characteristic angular 
dependence, despite possible interferences between them. The data demand very strongly the 
presence of an /o with mass 2105 ± 15 MeV. When it is dropped from the fit, gets worse by 
a very large amount, namely 4030. The description of the integrated 7777 cross section is then 
much worse, and there are also discrepancies with 7777 differential cross sections, illustrated with 
examples in Figs. 9(p)-(r). This resonance has the best determined mass and width of all the 
resonances observed in the present analysis, except for /4(2020). 

The mass agrees well with that observed by the E760 group [27] as a strong peak in rjrj in the 
reaction pp — > 777777° at two beam momenta of 3.1 and 3.5 GeV/c. The width we fit here, namely 
200 ± 25 MeV, is close to that fitted by E760: 203 ± 10 MeV. These parameters also agree well 
with the mass fitted to a 0+ peak in the 4tt channel in J/^ — )■ 7(47r) [28]. This resonance has 
also been identified as having J = in an analysis of Crystal Barrel data on pp — t/t/tt^ [29], 
where it appears as a strong peak in 7777; there, the determination of mass and width were not so 
precise, since the resonance appears at the top end of the available phase space. It is presently 
incorrectly listed by the Particle Data Group under /2(2150). 

We find a strong requirement for a further 0"'' resonance at 2320 MeV. If it is omitted, the 
fit gets visible worse, as is illustrated in Fig. 9(s); increases by 1257, a highly significant 
amount. 

There has been a report of a further 0+ resonance at 2020 MeV with a width of 400 MeV 
[30]. Adding this resonance to the fit, improves by 370, which we regard as significant. 
However, because it lies towards the bottom end of the available data for n^n^ and 7777 and 
because of interference with /o(2105), its mass and width have sizeable errors, M = 2005 ± 30 
MeV, F = 305 ± 50 MeV. 

4 Systematics of resonances masses 

On Fig. 10(a), we plot mass squared of 2+ resonances against their recurrence number. It is 
possible to construct an almost straight fine from /2(1270), /2(1565), /2(1910) and /2(2230).This 
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Figure 10: (a) Suggested trajectories for states (full line) and '^^2 (dashed) v. radial quantum 
number n; (b) likewise for 4+ states (full line), 3~ (dashed) and 1~ (dotted); (c) trajectories for 
3+ and 2~ states. Numerical values give masses in MeV. 



is illustrated by the full line on Fig. 10(a). We make use of one result as yet unpublished 
[31]. The mass of the /2(1565) has been obtained from an analysis of Crystal Barrel data on 
pp cucuTT^ at rest. From a fit with a Flatte form, the K-matrix mass is determined to be 
1598 ± n{stat) ± 9{syst) MeV. 

The /2(2020) and /2(2300) are readily placed on a parallel trajectory, shown by the dashed 
line of Fig. 10(a). The proximity of the /2(2020) to /4(2020) then suggests strongly that it is 
the qq n = 1 state and the /2(2300) is its radial excitation. The states on the full line of 
Fig. 10(a) are naturally interpreted as qq states. We interpret the shift in mass between the 
two lines as originating from the centrifugal barrier in the qq system. This provides an effective 
repulsion at small radii which shifts the resonance mass up more for L — 3 (^^2) than L — 1 

The slopes of the lines on Fig. 10(a) are 1.10 GeV^, with an error of ±0.03 GeV^. We find 
that it is possible to construct similar trajectories of identical slope for other quantum numbers. 
Fig. 10(b) shows corresponding trajectories for 4+, 3~ and 1~ states. There is a well defined 
shift in mass between p3(1960) and /4(2020), which may again be attributed to the effect of the 
centrifugal barrier, namely 50-60 MeV per unit of L. The dotted line on Fig. 10(c) shows a 
possible trajectory made up of p(1700) and p(2005). 

For = 2~, there is some indication of a discrepancy with a straight trajectory. However, 
there is evidence for an extra state 772 (1875) [32]. If all these 2~ states are confirmed, this is 
likely to be an intruder state. It has a mass roughly that expected for a 2~ hybrid [33]. 

From the analysis of rjn^n^ data, we can identify corresponding trajectories for 7 = qq 3+ 
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Figure 11: Regge trajectories for mesons with spins 1 to 6. Numerical values give masses in 
MeV. 



and 2~ states. They are shown on Fig. 10(c). For 3+, they are compatible with a trajectory 
of the same slope as other on Fig. 10, and with states approximately degenerate with 
and ^^2. The assignment of 0"*" resonances is presently controversial, so we do not show a 

trajectory for 0+. 

Fig. 11 shows Regge trajectories. The leading trajectory made up of p(770), /2(1270), 
P3(1690), /4(2020), p5(2295) and /6(2485) is well known. These resonances he close to a straight 
line versus s of slope 1.18 GcV^ per unit of spin. One can just discern a possible slight difference 
(~ 0.1 GcV^) between the trio 1^, 3" and 5~ and the even spin states. The slope is significantly 
greater than that of the lines on Fig. 10. Proof of this is the mass difference between /4(2020) 
and its grand-daughter /2(1910), likewise between p5(2295) and p3(2210) and the very clearly 
determined difference between /4(2020) and p3(1960). The centrifugal barrier in the qq system 
provides an explanation why high J states move up in mass. The intercept at zero mass with 
the vertical axis is at -0.67, slightly lower than the accepted value of -0.55 deduced from -nN 
scattering data [34]; however, the systematic shift between even and odd spins on this trajectory 
allows the possibility of moving the intercept up by 0.1. 

Fig. 11 also shows further parallel Regge trajectories. It is interesting that p(1700) and 
p(1450) depart significantly from straight-line trajectories through states of higher spin. Don- 
nachie and Clcgg [35] have argued that decay modes of p(1450) require it to be a hybrid or 
mixed with a hybrid. 



5 Flavour mixing 

We have explained in section 2 that there is a blurred distinction between (i) allowing phases 0j of 
decay channels to vary from or (ii) allowing flavour mixing through the angles Nonetheless, 
we feel it is worth summarising the fltted values of flavour mixing angles $ for 2"*" and 4"'" states. 
These are shown in Table 3. Errors on individual flavour mixing angles are typically ±5°. 
However differences in $ between different resonances are strongly correlated. These differences 
are what allow an acceptable fit to the data via small deviations from SU(3). Nonetheless some 
fiavour mixing is inevitable. If they are all set to zero, there is an unacceptable increase of 2646 
in x^- 

Table 3 also includes fitted values of corresponding ratios r of amplitudes for spins 1, 3,4 
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and 5. Errors are typically ±15%. One discerns a trend from positive values for low masses to 
negative values at high masses. As the mass rises, the momentum in the pp channel increases 
and the first zero of jL{kr) describing the pp channel moves to smaller radius. As it moves inside 
the radius where the resonance wave function peaks, one expects a change of sign of the ratios 
of amplitudes r between L = J ± 1, in qualitative agreement with the observations. However, 
because high mass resonances have several nodes in the wave function of the resonance, it is 
difficult to be quantitative for — 2+ and 1~ states at present. 



Resonance 




rj 


/2(1910) 


1.1 


2.89 


/2(2020) 


7.9 


1.62 


/2(2230) 


7.5 


-0.38 


/2(2300) 


-14.8 


-0.92 


/2(2620) 


(0) 


-0.83 


/4(2020) 


-26.1 


0.04 


/4(2300) 


22.9 


0.35 


P5(2295) 




(0) 


P5(2500) 




0.24 


P3(1960) 




0.04 


P3(2210) 




0.35 


P3(2300) 




2.0 


Pi(2005) 




5.0 


Pi(2165) 




-0.22 



Table 3: Values of flavour mixing angles $ and ratios of amplitudes rj — |/l=j+i|/|/l=j-i| as 
s — >■ oo. 

A possibility is that the flavour mixing arises purely from the overlap with neighbouring ss 
states. The ss partner of /2(1565) could be the 2"*" state reported by the LASS group at 1950 
MeV [36] decaying to K*K*, in which case mixing with /2(1910) would be natural. On the 
other hand, the LASS group may simply be observing the decays of /2(1910). The ss partner of 
/2(1910) itself is expected around 2150-2250 MeV. It is possible to identify the states observed 
in (j)(f) as the ss partners of /2(1910) and /2(2020). The state decaying to the 00 S-wave, 
reported by Etkin et al. [37], also peaks at about 2150 MeV; however, their K-matrix analysis 
assigns it a mass of 2020 McV, because of a strong threshold effect. The JETSET group have 
likewise reported a signal in pp — )> (pep peaking at ~ 2180 MeV, possibly the same object [38]. 
Prokoshin has reported a peak in rjrj at 2175 ± 20 MeV [39], again possibly the same resonance. 
Recent data of the Omega group on central production of K'^K~ reveals a peak at 2150 MeV 
[40]. 

Etkin et al. have also reported a peak in the 00 D-wave at 2300-2340 MeV. This makes a 
natural candidate for the ss partner of /2(2020). 

Production of the 00 resonances from initial pp and tttt states requires explanation. It could 
arise if resonances arc strongly mixed between qq and ss. However, the small flavour mixing 
angles we observe do not point that way. 
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A second clear possibility is that the mixing is with the 2+ glueball predicted in this mass 
range, or slightly above, by Lattice QCD calculations. There is evidence from the Crystal 
Barrel data on pp — )■ r]r]7r^ for a broad 2+ state in rir] with a mass of 1980 ± 50 MeV and a width 
r = 500 ± 100 MeV [29]. There is similar evidence from two other sources. Firstly, a broad 
2+ contribution to the Air channel is observed in central production data at small pt [30]. Its 
mass is 1920 ± 20 MeV with a width of 450 ± 60 MeV. Secondly, there are data from BES on 
J/^ — > jAn [41] and J/^ — )■ '~fK*K* [42], which both require a broad 2"^ signal peaking at about 
2000 MeV. The presence of this wide 2+ signal in J/^ radiative decays and central production 
data is suggestive of mixing between a glueball and neighbouring qq states to produce a broad 
state. 

We have tried inserting a broad /2(1980) with F = 500 MeV into our analysis, in addition to 
the four states discussed above. There is a modest improvement in of ~ 100, but little 
visible change to the quality of the fits. This is not sufficient to confirm the possible presence of 
a broad 2+ background in the present data, because of the complexity of possible interferences 
with /2(1910) to /2(2300). 

5.1 /o(2105) 

One very striking result of the present analysis is that /o(2105) is far from being a normal 
qq state. Its flavour mixing angle $ is +(59 — 64)° if /o(2005) is excluded from the fit, or 
+ (68 — 71.6)° with it included. That is, it behaves predominantly as an ss state. Allowing 
for interferences with /o(2005), it makes up (4.6 ± 1.5)% of the tt^tt^ = 0^ intensity and 
(38 ±5)% of that for rji]. The branching ratio to rji]' is not well determined, because of the low 
statistics in this channel. Our best estimate of amplitude ratios is as follows: 

ttV : 7777 : 7777' = 0.71 ± 0.17 : 1 : -0.85 ± 0.45. (15) 

For an unmixed qq state, the ratio expected between n^n^ and 7777 is 0.8~^, i.e. 2.44. The strong 
production of an ss state in pp annihilation is clearly anomalous. How can this be explained? 

It could be a second glueball. A glueball would have a mixing angle of +37°. The fiavour 
mixing angle we observe then required some mixing with a nearby ss state; for ideal mixing, 
ss states are members of SU(3) singlets, so such mixing is plausible. The latest Lattice QCD 
calculations [43] predict a second 0+ glueball with a mass ratio to the first of 1.54 ± 0.11. If we 
assign /o(1500) and /o(2105) as the two glueballs, the mass ratio is 1.40. 

The alternative scenario, developed by Anisovich and Sarantsev [44], is that the glueball 
mixes strongly with neighbouring qq and ss states to generate a broad state which contains 
~ 50% of the glueball in its wave function. Such a broad state may well overlap the mass range 
of the /o(2105). The /o(1500), /o(1770) and /o(2105) are all seen strongly as peaks in J/* 
radiative decay to An [28]. This suggests that /o(2105) may be a mixed state made by mixing 
of the glueball with a nearby ss state. It is important to look for it in decays to KK. 

6 Conclusions 

A combined analysis of the Tr^TT"*" data of Eisenhandler et al. and Hasan et al. with Crystal 
Barrel data for 7r°7r°, 7777 and 7777' leads to a secure partial wave analysis and the identification of 
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most of the expected qq resonances in this mass range, except for two 1~ states. The polarisation 
data play a vital role in separating states with L = J ± 1. The resulting fit is rather close to 
that obtained previously by Hasan and Bugg [9]. However, it is now possible to identify clearly 
four 2+ resonances. These have small flavour mixing and are probably to be identifled with the 
qq and states expected in this mass range. 

The systematics of resonances on Fig. 10 establishes a valuable guide to masses to be expected 
to qq states. This should help in identifying intruder states of exotic character. 

The /o(2105) make a large contribution to the rjrj and rjrj' data. It is not to be identifled with 
a simple qq state, because of its +(59-71.6)° flavour mixing angle. This requires it to have a 
large ss and/or exotic component. 
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